Time-resolved measurement of Landau-Zener tunneling in periodic potentials 
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We report time-resolved measurements of Landau-Zener tunneling of Bose-Einstein condensates in 
accelerated optical lattices, clearly resolving the step-like time dependence of the band populations. 
Using different experimental protocols we were able to measure the tunneling probability both in the 
adiabatic and in the diabatic bases of the system. We also experimentally determine the contribution 
of the momentum width of the Bose condensates to the width of the tunneling steps and discuss 
the implications for measuring the jump time in the Landau-Zener problem. 

PACS numbers: 03.65.Xp, 03.75.Lm 
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Tunneling is one of the most striking manifestations of 
quantum behaviour and has been the subject of intense 
research both in fundamental and applied physics 
While tunneling probabilities can be calculated accurately 
even for complex quantum systems and have an intuitive 
interpretation as statistical mean values of experimental 
outcomes, the concept of tunneling time and its compu- 
tation are still the subject of debate even for simple sys- 
tems Q . How much time a quantum system spends in a 
classically forbidden area, e.g., inside a barrier separating 
two potential wells, has been theoretically calculated for 
a wide range of physical systems [3, 0] and measured in 
recent experiments [f|. A well-known example of a sim- 
ple tunneling problem is Landau-Zener (LZ) tunneling, 
in which a quantum system tunnels across an energy gap 
at an avoided crossing of the system's energy levels. In 
recent years, LZ tunneling has been used as a paradigm 
for studies of the the quantum Zeno effect (controlling 
the decay by repeated measurements of the systems 0]), 
and more generally in the observation of deviations from 
a purely exponential decay of a quantum system consist- 
ing of just a few energy levels. An example of the latter 
was investigated by Raizen and co-workers where 
the decay in a Wannier-Stark system was shown to de- 
viate from an exponential law at short times. Moreover, 
LZ tunneling provides a building block for the quantum 
control of complex many-body systems [ll[ . 

According to the adiabatic theorem an infinitely slow 
change in the Hamiltonian of a quantum system leads 
to an adiabatic following of the instantaneous eigenstate, 
whereas a controlled time dependence can cause transi- 
tions to other eigenstates. In the case of a linear time de- 
pendence in a Hamiltonian describing two crossing eigen- 
states that are repelled by a coupling energy, the prob- 
ability of a transition between these eigenstates is de- 
scribed by the celebrated Landau-Zener model [l2l [T^]. 
While in most treatments of LZ-tunneling only the 
(asymptotic) tunneling probability is considered, the full 
time dependence of the LZ dynamics as determined by 
the time evolution of the quantum-mechanical wavefunc- 
tion also plays an important role, especially if the asymp- 



totic limit assumed in LZ theory is not reached in a sin- 
gle transition. The characteristic time required for the 
tunneling between initial and final quantum states to 
be completed 0, Q d is sometimes referred to as the 
"jump time" . 

In this Letter, we present experimental results, backed 
up by numerical simulations, on the LZ dynamics and 
on the jump time for a Wannier-Stark system realized 
with ultracold atoms forming a Bose-Einstein condensate 
(BEC) inside an optical lattice [la, ll7| • This system gives 
us an optimal control of single as well as sequences [18| of 
LZ processes. In contrast to the above mentioned experi- 
ments with cold atoms 0, H| our BEC has an initial width 
in momentum space that is much smaller than the char- 
acteristic momentum scale of the problem given by the 
momentum width ps = 2p rec = 2-K%/di, of the first Bril- 
louin zone of a periodic potential with lattice constant 
dr. This possibility of engineering the initial state for the 
LZ process in momentum space enables us to observe the 
time dependent dynamics with a non-exponential decay 
in the survival probability for single or multiple LZ cross- 
ings as well as oscillations in the transition probabil- 
ity after the system has passed the avoided crossing (as 
predicted for LZ tunneling in atomic Rydberg states [2l| 
and experimentally observed in a wave-optical two-level 
system [22|), the only limitation being the initial mo- 
mentum width of the condensates and nonlinear effects. 
This possibility of engineering the initial state for the LZ 
process in momentum space enables us to observe the 
time dependent dynamics with a non-exponential decay 
in the survival probability for single or multiple LZ cross- 
ings as well as oscillations in the transition probabil- 
ity after the system has passed the avoided crossing, the 
only limitation being the initial momentum width of the 
condensates and nonlinear effects. Such oscillations were 
predicted for LZ tunneling in atomic Rydberg states [2lJ 
and experimentally observed in a wave-optical two-level 
system [22j]). Our experiments are similar to recent stud- 
ies of LZ transitions at avoided crossings in the energy 
levels of a solid-state artificial atom [23|], but the high 
level of control over the light-induced periodic potential 
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FIG. 1: Time-resolved measurement of LZ tunneling, (a) 
Experimental protocol (shown in the band-structure repre- 
sentation of energy E(q) versus lattice quasimomentum q). 
Left: The lattice is accelerated, (partial) tunneling occurs. 
Right: The acceleration is then suddenly reduced and the lat- 
tice depth increased so as to 'freeze' the instantaneous popu- 
lations in the lowest two bands; finally, further acceleration is 
used to separate these populations in momentum space, (b) 
Experimental results for Vb = 1 E rcc and _Fo = 0.383, giving 
Tb = 0.826 ms. The solid and dashed lines are a numerical 
simulation of our experimental protocol and an exponential 
decay curve for our system's parameters, respectively. 

also allowed us to measure the tunneling dynamics in 
different eigenbases of the system's Hamiltonian. 

In our experiments we created Bose-Einstein conden- 
sates of 5 x 10 4 rubidium-87 atoms inside an optical 
dipole trap (mean trap frequency around 80 Hz). A 
one-dimensional optical lattice created by two counter- 
propagating, linearly polarized gaussian beams was then 
superposed on the BEC by ramping up the power in the 
lattice beams in 100 ms. The wavelength of the lattice 
beams was A = 842 nm, leading to a sinusoidal poten- 
tial with lattice constant cIl = A/2 = 421 nm. A small 
frequency offset between the two beams could be intro- 
duced through the acousto-optic modulators in the setup, 
allowing us to accelerate the lattice in a controlled fash- 
ion. 

The time-resolved measurement of LZ tunneling was 
then effected as follows [see Fig. 1 (a)]. After loading 
the BEC into the ground state energy band of an optical 
lattice of depth Vb as described above, the lattice was ac- 
celerated with acceleration auz for a time £lz by chirping 
the frequency offset between the lattice beams, resulting 



in a corresponding force -FLz on the atoms in the rest 
frame of the lattice. The lattice thus acquired a final ve- 
locity v = olz^lz- During t^z the quasimomentum of the 
BEC swept the Brillouin zone, and at multiples of half 
the Bloch time T B = 2irh(Ma LZ d L y 1 (where M is the 
atomic mass) i.e. at times t = (n+l/2)Ts (n = 0, 1, 2, ...) 
when the system was close to the Brillouin zone edge, 
tunneling to the upper band became increasingly likely. 
At time t — t^z the acceleration was abruptly reduced to 
a smaller value a scp and the lattice depth was increased 
to V scp in a time t ramp <C Tb. These values were cho- 
sen in such a way that at time t = t\jz the probability 
for Landau-Zcncr tunneling from the lowest to the first 
excited energy band dropped from between « 0.1 — 0.9 
(depending on the initial parameters chosen) to less than 
w 0.01, while the tunneling probability from the first ex- 
cited to the second excited band remained high at about 
0.95. This meant that at t = £lz the tunneling process 
was effectively interrupted and for t > t^z the measured 
survival probability P(t) — No/N tot (calculated from the 
number of atoms N in the lowest band and the total 
number of atoms in the condensate AT to t) reflected the 
instantaneous value P(t = 4lz)- 

The lattice was then further accelerated for a time t sep 
such that a S episcp ~ 2np rcc /M (where typically n = 2 
or 3). In this way, atoms in the lowest band were ac- 
celerated to a final velocity v s=s 2np roc /M, while atoms 
that had undergone tunneling to the first excited band 
before t = £lz underwent further tunneling to higher 
bands with a probability > 0.95 and were, therefore, no 
longer accelerated. At time t scp the lattice and dipole 
trap beams were suddenly switched off and the expanded 
atomic cloud was imaged after 23 ms. In these time-of- 
flight images the two velocity classes and 2np TCC /M 
were well separated and the atom numbers Nq and JV"t t 
could be measured directly. Since the populations were 
"frozen" inside the energy bands of the lattice, which 
represent the adiabatic cigenstates of the total Hamilto- 
nian of the system, this experiment effectively measured 
the time dependence of the LZ survival probability P a in 
the adiabatic basis. The result of a typical measurement 
is shown in Fig. 1 (b). One clearly sees two "steps" at 
times t = 0.5 Tb and t = 1.5 Tb, which correspond to 
the instants at which the atoms cross the Brillouin zone 
edges, where the lowest and first excited energy bands ex- 
hibit avoided crossings. For comparison, the result of a 
numerical simulation (integrating the linear Schrodingcr 
equation for the experimental protocol) as well as an ex- 
ponential decay as predicted by LZ theory are also shown. 

The LZ tunneling probability can be calculated by con- 
sidering a two-level system with the adiabatic Hamilto- 
nian 

AE . . 

H a = H d + V = atcTy + —<?x, (1) 

where <7j are the Pauli matrices. The eigenstates of the 
diabatic Hamiltonian H^, whose eigenenergies vary lin- 
early in time, are mixed by the potential V characterized 
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FIG. 2: (a) LZ survival probability in the adiabatic basis for 
a fixed force Fo — 1.197 and different lattice depths (filled 
squares: Vo = 2.3E rcc ; open circles: Vq = 1.8E rec ; open 
squares: Vo = lE rcc ; filled circles: Vo = 0.6E rcc ). The dashed 
lines are sigmoid fits to the experimental data. Inset: Survival 
probability in both the adiabatic (open squares) and diabatic 
(filled triangles) bases for Vo = 1 E roc and Fo = 1.197. (b) 
Step height ft as a function of the inverse adiabaticity pa- 
rameter I/7 for varying lattice depth and Fo = 1.197 (open 
symbols), and for varying force with fixed Vo = 1.8 E rcc (filled 
symbols) . The dashed line is the prediction of Eq. (3) for the 
LZ tunneling probability. 



by the energy gap AE. Applying the Zener model [l3[ to 
our case of a BEC crossing the Brillouin zone edge leads 
to a band gap AE = Vq/2 and to a = 2v lec Mai,z = 
2F E? ec /(nh), with E lec = H 2 tt 2 /(2Md 2 ) the recoil en- 
ergy and Fq — Molz^l/ E rec the dimensionless force. 
The limiting value of the adiabatic and diabatic Lanadau- 
Zener survival probabilities (for t going from — 00 to +00) 
in the eigenstates of H a and Hd, respectively, is 



P a (t -> +00) = 1 - P d (t - 
where 

Plz = 



-00) = 1 - F\ 



LZ, 



(2) 



(3) 

4ha(AE)- 2 is the 



with the adiabaticity parameter 7 
standard LZ tunneling probability [24l |. 

In order to test whether this tunneling probability cor- 
rectly predicts the height of the step corresponding to a 



single LZ tunneling event, we performed the experiment 
described above for a variety of values of V and F. Fig- 
ure 2(a) shows the first tunneling step for different lattice 
depths Vo, measured in units of F rec at a given LZ accel- 
eration determined by Fo- The steps can be well fitted 
with a sigmoid function 



Pait) = 1 - 



1 + exp((td - tJ/AtLz) ■ 



(4) 



where to is the position of the step (which can deviate 
slightly from the expected value of 0.5 Tb, e.g. due to a 
non-zero initial momentum of the condensate), h is the 
step height given by Eq. ([3]) and At^z represents the 
width of the step. As expected, the step height h follows 
the LZ result for the tunneling probability (see Fig. 2(b)). 

We note here that while the experimental protocol 
described above measures the LZ tunneling probabil- 
ity in the adiabatic basis, it is also possible to make 
analogous measurements in the diabatic basis of the un- 
perturbed free-particle wavefunctions (plane waves with 
a quadratic energy-momentum dispersion relation) by 
abruptly switching off the lattice and the dipole trap after 
the first acceleration step. In this case, after a time-of- 
flight the number of atoms in the v = and v = 2p Iec /M 
momentum classes are measured and from these the sur- 
vival probability (corresponding to the atoms remaining 
in the v — velocity class relative to the total atom 
number) is calculated. The inset of Fig. 2(a) (filled tri- 
angles) shows such a measurement. Again, a step around 
t = 0.5 Tb is clearly seen, as well as strong oscillations 
for t > 0.5 7s. While (weaker) oscillations can also 
be measured in the adiabatic basis (see the results for 
Vo = 2.3E rcc in Fig. 2(a)), they are much stronger and 
visible for a wider range of parameters in the diabatic 
basis, as expected from theoretical calculations [Hj]. 

The widths A^lz corresponding to the steps shown in 
Fig. 2(a) should, according to our interpretation, re- 
flect the "tunneling time" or "jump time" for LZ tunnel- 
ing AiLz = At>Lz/ftLZ during which the probability of 
finding the atoms in the lowest energy band goes from 
P a (t = 0) = 1 to its asymptotic LZ value Plz- In the 
theoretical literature this tunne ling time has been calcu- 
lated by several authors 0, [T3, Il5l | Vitanov [l5[ defines 
the jump time in the adiabatic basis as 



jump 



Pg(t= +Q0) 



(5) 



where P' a (t = to) denotes the time derivative of the 
tunneling probability P a (t) evaluated at the crossing 
point of H a . A sigmoidal function for P a (t) leads to 
T iump _ 4A£ LZ> p or l ar g e values of 7, which is the regime 
of our experiments, the theoretical jump time is given by 
rl ump /T B w (Vo/F rcc )/8. From our sigmoidal fits we 
retrieve T^ Lmp /T-q w 0.15 — 0.35 (corresponding to abso- 
lute jump times between 50 /is and 200 /is), whereas the 
theoretical values for our experimental parameters are in 
the region of 0.1 — 0.15. We interpret this discrepancy 
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as being due the fact that in our experiment the con- 
densate does not occupy one single quasimomentum but 
is represented by a momentum distribution of a finite 
width Ap/ps ^ 0.1 due to the finite number of lattice 
sites (around 50) it occupies and the effects of atom-atom 
interactions. Also, during the acceleration process dy- 
namical instabilities can further broaden the momentum 
distribution. 

In order to test the dependence of the measured step 
width on the width of the momentum distribution of 
the condensate, we created initial distributions of differ- 
ent widths using a dynamical instability in a controlled 
way [26|. The condensate was loaded into a lattice mov- 
ing at a finite velocity corresponding to quasimomentum 
q = —O.Sps and held there for a time up to 3 ms. Dur- 
ing this time the dynamical instability associated with 
the negative effective mass at q = — O.Spb led to an in- 
crease in the momentum width of the condensate. After 
this preparatory stage, the LZ dynamics was measured 
as described above and the width of the tunneling step 
was measured [see Fig. 3(a)]. As expected, the larger the 
initial momentum width of the condensate, the larger the 
step width [Fig. 3(b)]. This was also confirmed by a nu- 
merical integration of the Schrodinger equation in which 
the momentum width of the condensate was varied by 
changing the initial trap frequency. The simulation also 
showed that for a vanishing momentum width, the step 
width still remains finite and in that limit directly reflects 
the jump time given by Eq. ([5]). 

In summary, we have measured the full dynamics of 
a single LZ transition of matter waves in an accelerated 
optical lattice. In both the adiabatic and diabatic bases 
the step-like behaviour as well as oscillations of the sur- 
vival probability were clearly seen and agree with the- 
oretical predictions. While there is quantitative agree- 
ment between the measured and theoretically expected 
step heights, the jump times calculated from the step 
widths are larger by ~ 50 - 100% than the theoretical 
results because of the finite initial momentum widths of 
the condensates in our experiments. 

In future investigations one could reduce the initial mo- 
mentum width by using, e.g., appropriate trap geometries 
or by controlling the nonlincarity through Fcshbach reso- 
nances. This would allow one to obtain a more quantita- 
tive measurement of the jump time and enable a compari- 
son with theoretical results related to the minimum time 
for a single LZ crossing limited by fundamental quan- 
tum (or wave, see (22[) mechanical properties [27|. Also, 
clearer observations of the short-time oscillations (hap- 



pening on even shorter time scales) , whose signatures can 
be seen in Fig. 2(a) should be possible in this way. Our 
method can also be used to study multiple LZ crossings, 
e.g., in order to observe Stuckelberg oscillations. 
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FIG. 3: LZ transition for different momentum widths of the 
condensate, (a) Survival probability for Ap/2p rcc — 0.2 (filled 
squares) and Ap/2p rcc = 0.6 (open squares). The solid and 
dashed lines are the results of a numerical simulation and of 
a sigmoid fit, respectively, (b) Step width AtLz of the LZ 
transition as a function of the momentum width of the con- 
densate. The solid line is the result of a numerical simulation. 



129/1) and the Global Networks Mobility Measures. We 
thank Martin Holthaus and Tilman Esslinger for discus- 
sions and comments on the manuscript. 



[1] M. Razavy, Quantum Theory of Tunneling, World Scien- 
tific, Singapore (2003). 

[2] L.S. Schulman, Lect. Notes Phys. 734, 107 (2008). 

[3] E.H. Hauge and J.A. St0vneng, Rev. Mod. Phys. 61, 917 
(1989). 

[4] M. Biittiker, R. Landauer, and Th. Martin, Rev. Mod. 



Phys. 66, 217 (1994). 
[5] P. Eckle et al., Science 322, 1525 (2008). 
[6] P. Facchi and S. Pascazio, J. Phys. A 41, 493001 (2008). 
[7] SR. Wilkinson et al., Nature (London) 387 575 (1997). 
[8] Q. Niu and M.G. Raizen, Phys. Rev. Lett. 80 3491 

(1998). 



5 



[9] M. C. Fischer, B. Gutirrez-Medina, and M. G. Raizen, 
Phys. Rev. Lett 7, 040402 (1999). 

[10] J.M. Raimond, M. Brune, and S. Haroche, Rev. Mod. 
Phys. 73, 565 (2001). 

[11] G.E. Santoro, R. Martonak, E. Tosatti, and R. Car, Sci- 
ence 295, 2427 (2002). 

[12] L.D. Landau, Phys. Z. Sowjetunion, 2, 46 (1932). 

[13] C. Zener, Proc. R. Soc. London, Ser. A 137, 696 (1932). 

[14] K. Mullen, E. Ben- Jacob, Y. Gefen, and Z. Schuss, Phys. 
Rev. Lett. 62, 2543 (1989). 

[15] N. V. Vitanov, Phys. Rev. A 59, 988 (1999). 

[16] B.P. Anderson and M.A. Kasevich, Science 27 282 
(1998). 

[17] O. Morsch, J.H. Miiller, M. Cristiani, D. Ciampini, and 

E. Arimondo, Phys. Rev. Lett. 87, 140402 (2001). 
[18] C. Sias et al., Phys. Rev. Lett. 98 120403 (2007). 



[19] Ref. [20I ] studied the influence of atomic interactions in 
BEC samples on the non-exponential decay in LZ se- 
quences. 

[20] S. Wimberger et al, Phys. Rev. A 72 063610 (2005). 
[21] J.R. Rubmark, M.M. Kash, M.G. Littman, and D. 

Kleppner, Phys. Rev. A 23, 3107 (1981). 
[22] D. Bouwmeester et al, Phys. Rev. A 51 646 (1995). 
[23] D.M. Berns et al., Nature (London) 455, 51 (2008). 
[24] Ref. [25| derived a more refined Zener model for ultracold 

atoms in an optical lattice. 
[25] M. Holthaus, J. Opt. B:Quantum Semiclass Opt. 2, 589 

(2000). 

[26] M. Cristiani et al, Opt. Express 12, 4 (2004). 
[27] Y. Makhlin, and P. Hako nen, Phys. Rev. Le tt. 96, 187002 
(2006); T. Caneva et al. larXiv:0902.4193l 



